Current driven defect unbinding transition in an XY ferromagnet 
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A Keldysh-contour effective field theory is derived for magnetic vortices in the presence of current 
flow. The effect of adiabatic and non-adiabatic spin transfer torques on vortex motion is highlighted. 
Similarities to and differences from the superconducting case are presented and explained. Current 
flow across a magnetically ordered state is shown to lead to a defect-unbinding phase transition 
which is intrinsically nonequilibrium in the sense of not being driven by a variation in effective 
temperature. The dependence of the density of vortices on the current density is determined. 
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INTRODUCTION 



The nonequilibrium physics of strongly interacting 
quantum systems is a topic of fundamental theoreti- 
cal importance. A question of particular current inter- 
est is the behavior of quantum critical systems under 
nonequilibrium conditions. Examples include phase tran- 
sitions induced at temperature T=0 in the presence of 
a nonequilibrium drive*— as well as quantum coarsening 
and quench problems of a system prepared or maintained 
in a far- from equilibrium state In many of the cases 
of nonequilibrium quantum criticality studied to date, 
the dominant physics is that the nonequilibrium drive 
acts to produce a noise term (typically delta-correlated 
at long scales) in the equation of motion for the criti- 
cal fluctuations i This noise term increases fluctuations 
and destroys order, similarly to temperature. Differences 
from equilibrium physics of course appear, for example 
as violations of the fluctuation-dissipation theorerr&H or 
as subleading corrections to scaling^^ Other situations 
can occur too. For example, Dalle Torre et al showed that 
in some circumstances departure from equilibrium can 
lead to a noise with long ranged correlations, which leads 
to larger differences from equilibrium physics.— In some 
models, excitations propagate ballistically so the concept 
of effective temperature does not apply and nonequilib- 
rium phase transitions arise for different reasons 

In this paper we investigate an alternative route to 
nonequilibrium quantum criticality. We show that in 
continuous-symmetry magnets an applied current can 
lead to a current-driven defect unbinding. This has no 
equilibrium analogue and therefore constitutes an intrin- 
sically nonequilibrium phase transition, quantum in the 
sense that it is driven at T—0. Of course, even at T=0 
a driving field will introduce a noise, but analogously 
to the nonequilibrium classical phase transitions occur- 
ring e.g. in sheared liquids^ the noise is not the cru- 
cial parameter. The mechanism we discuss is generic, 
applying (with some difference of details) for any combi- 
nation of dimensionality and order parameter symmetry 
for which topological defects can be denned. For definitc- 
ness we present specific results in the case of a quasi two- 



dimensional XY magnet, for which the defects involved 
are vortices. In this respect our work complements the 
treatment of current-driven quantum criticality in our 
previous papers^^ which treated Ising and Heisenberg 
symmetries and found that the transition was driven by 
a current-induced noise. 

That an applied current can lead to a force on a mag- 
netic defect was previously known.— ~— The transition we 
find is analogous to defect-unbinding transitions driven 
by supercurrent flow in superfluidsi 19 ' 20 The new contri- 
butions of this paper include a derivation from microscop- 
ies which is somewhat different from (although basically 
in agreement with) previous workj2!~— a discussion of 
the mapping between the superconducting and magnetic 
situations, some additional insights into the relative im- 
portance of the different mechanisms by which Galilean 
invariance is broken and explicit results for the defect 
creation rates. 

Our approach is to write a general interacting fermion 
model in the presence of a nonequilibrium drive as a 
Keldysh-contour path integral, introduce magnetic de- 
grees of freedom via a Hubbard-Stratonovich transfor- 
mation of the fermion-fermion interactions and by inte- 
grating out the fermions obtain a nonequilibrium spin 
model. We use a method presented by Schula 2 ^ to intro- 
duce rotational degrees of freedom and use an expansion 
in gradients and quantum fluctuations to obtain an action 
in terms of rotational spin excitations. The mean field ap- 
proximation to this action is the Landau-Lifshitz-Gilbert 
equations. By representing magnetic configurations in 
terms of topological (vortex) coordinates^ we obtain an 
action for the topological excitations, which is then ana- 
lyzed via quantum critical and instanton methods. 

The rest of this paper is organized as follows. Sec- 
tion [IT] outlines the approach we use, Sections IIIII and ITVl 
present a derivation of the basic dynamical effects, Sec- 
tion [V] summarizes the resulting action in spin coordi- 
nates, Section IVII derives the action for topological de- 
fects, Section rVIII describes the current-driven vortex un- 
binding transition and Section I Villi is a summary and 
conclusion. 
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II. FORMALISM 
A. Overview 

We are interested in the effect of current drive on the 
nucleation and motion of topological defects in the mag- 
netization field of a metallic ferromagnet. Although our 
formalism can be applied to the case of any topologi- 
cal defect, we present explict results for vortices: the 
topological defects in a two dimensional easy-plane mag- 
net. Because topological defects are large scale structures 
in the spin field, they obey equations of motion which 
may be deduced from semiclassical arguments based 
on the Landau-Lifshitz-Gilbert equation j 16 i 17 i 22 Denot- 
ing the position of defect i by Xi and the sign of its 
vorticity by r\ = ±1, we may represent its equation of 
motion as 

tjCJozxX.+IdVqX, = r 1 CzxJ+I D V+^F ij +Ci (1) 

3 

with Fij the force due to the intervortex interaction 
(including both spin- wave- mediated and dipolar terms). 
The terms J and V arise due to current flow, while 
and T>o representing the effects of spin conserving and 
non-conserving processes on the temporal dynamics. Id 
and C are constants and Q is a fluctuating field which on 
scales longer than the mean free path may be taken to 
be delta correlated, with 

«?(*)#(*') = 2I D N XX T* SijS^Sft - (2) 

Here N xx is related to the many-body density of states 
and T* is an effective temperature arising from the de- 
partures from equilibrium. 

In this paper we employ a Keldysh-contour based 
nonequilibrium approach to calculate the coefficients in 
Eq. in terms of the underlying properties of the elec- 
trons (in particular the transport and spin relaxation 
times) and the applied current, determining the rela- 
tive magnitudes of the different terms and explicating the 
similarities and differences between topological defects in 
superconductors and in magnets. 

The formalism is based on a phase angle representa- 
tion convenient for carrying out a gradient expansion. It 
enables a quantitative understanding of the origin of the 
breaking of Galilean invariance, of the roles of spin con- 
serving and non-conserving scattering mechanisms, and 
of the variations of parameters as the magnetic critical 
point is approached. We show that in the generic sit- 
uation when the spin non-conserving scattering is small 
compared to the spin conserving scattering, then the vor- 
tex Hall angle is typically near to 90° (in contrast to the 
result for superconductors), but the drag contribution 
to the force on a vortex is much larger than the trans- 
verse ('Magnus') force contribution (also different from 
the usual case of vortices in superconductors). These re- 
sults, in combination with the effects of the long-ranged 



dipole interaction, lead to quantitative differences in the 
current-induced nucleation of vortices. 



B. Basic Assumptions 

We consider a system of interacting fermions and make 
the Fermi liquid assumption that the dynamical effects of 
the electron-electron interaction may be neglected at the 
scales of interest here. The system is then characterized 
by a band structure, which we represent via a hopping 
matrix ti-j connecting unit cells i and j, and a local 
magnetization 

fh(r,t) — h(r,t)m(f,t) (3) 

of direction n and magnitude m. 

The physics also involves a self energy £ which en- 
codes the effects of scattering processes. We will distin- 
guish two sources of scattering: spin conserving processes 
corresponding e.g. to scattering off of spin-independent 
random potentials and spin non-conserving processes: 

Stot = S cons + Xat-c* (4) 

We will assume that the spin non-conserving processes 
arise from electron exchange with a spin-unpolarized 
reservoir which does not change when a current is ap- 
plied to the systemJ^. Thus Sjv-c is proportional to 
the unit matrix in spin space and is independent of ap- 
plied current. Because for example the density of states 
may be spin dependent S cons may have structure in spin 
space. For simplicity we will assume that the magnetiza- 
tion varies sufficiently slowly in space that the electrons 
always relax to the local magnetization; in this case S cons 
is the sum of a term proportional to the unit matrix and 
a term proportional to m • a with, in particular, no terms 
proportional to gradients of m. Where an explicit model 
is necessary we will take E cons to arise from randomly 
placed point-like impurities. 

We suppose that the system is driven out of equilib- 
rium by an externally imposed current which is weak 
enough to be described by linear response theory. Since 
we are interested in length scales long compared to the 
spin relaxation length, any spin polarized current in- 
jected at the boundary will quickly relax. The relevant 
currents are thus induced by applying an electric field, 
E. In the presence of a non-vanishing magnetization the 
resistivity p will be a matrix in spin space so the current 
density J = p~ l E will have both a charge and a spin 
component. Our assumption that electrons are relaxed 
to the local magnetization direction means that the spin 
direction is aligned to the local magnetization so that 
in the basis aligned with the local magnetization p has 
components p a which may differ for up and down spins 
and (denoting the Cartesian component of a vector by a 
roman symbol) 

J a = J Q 1 + J a M h-S (5) 
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In the Keldysh two contour Greens function method 
the current flowing between sites i and j is determined 
by the Keldysh component G K of the electron Greens 
function via 



J^=-iG K (i,j;t,t)J^ 



with current operator 



J'l—j l&i 



-j (n 



(6) 



(7) 



Here (r*j — fj) is the vector connecting sites i and j and 
in momentum space J — dek/dk = Vk 

The retarded and advanced components of the Greens 
function and self energy change only to second order 
in departures from equilibrium, so within the linear re- 
sponse approximation, departures from equilibrium may 
be described by nonequilibrium terms in the Keldysh 
components of the electron Greens function G and self 
energy S. Thus we have e.g. G A ' -> Gf q + G% eq , 
with G^eq calculated to linear order in the applied field. 
The nonequilibrium terms are discussed in detail in sec- 
tion IIV1 here we note simply that in a simple electron- 
ically three dimensional model with k 2 /2m dispersion, 
magnetization aligned along the spin-z axis and rela- 
tively weak scattering, the formalism leads to the familiar 
Drude formula 



Pa 



(8) 



with r^. ot the possibly spin-dependent scattering rate ob- 
tained from 7mS fot (cj = 0) and n a the total number of 
electrons with spin aligned (or anti-aligned) to the local 
magnetic field. 

Also, as discussed in Ref. H3 an applied electric field 
leads to a noise term corresponding to an effective tem- 
perature 



T. 



eff 



eEl* 



(9) 



related to the voltage dropped over an appropriate length 
scale /*. The differences in physics between the model 
considered here and the model studied in Ref. [HI will 
lead to a different expression for /*. 



C. Nonequilibrium terms in Greens function 

In this subsection we sketch a derivation of the 
nonequilibrium terms in the electron Greens function 
and self energy. The treatment follows Ref. 13. Not- 
ing that the retarded and advance self energies are un- 
changed to linear order in applied fields, we parametrize 
the change in the Keldysh components of G, X via a dis- 
tribution function h related to the Greens function via 
G K = G R h - hG A and similarly for S. We then follow 
the usual steps to obtain a quantum kinetic equation for 
h. In writing the kinetic equation we make use of our 



assumption that the electron distribution relaxes to the 
local magnetization so that we have one kinetic equation 
for each spin direction. We further observe that in the 
model considered here, impurity scattering may relax the 
momentum dependence of the distribution function but 
not the energy dependence. We write h = h eq + h A + h s 
with h eq (x) = 1 — 20(— x) and nonequilibrium parts h A 
odd and hs even under k O —k. Separating the even and 
odd parts of the kinetic equation we obtain (not denoting 
the spin direction explicitly) 



^ (E ■ v) + 2Y tot h A = -25{x) (E ■ v) 



and 



dh A 
dx 



E-V + 2T 



N-cn-s 



(10) 



(11) 



Solving and restoring spin indices explicitly we find 



w = - 



and 



E-v VT 



N-CL tot 



per 
1 tot 



E ■ V 



Exp 



2\x\y/T 



N-C^tot 



E ■ v 



h% = -2 



sgnix) 



Exp 



2 \x\ y/T N - C T° t 



(12) 



(13) 



These expressions differ from those obtained in Ref. LL3| 
because in that work T tot = Tjv_c- 

For small E the term in parentheses in Eq. (|12p is 
just an integral representation of the delta function 5{x). 
To an accuracy sufficient for the leading terms in the 
semiclassical expansion fSec. IIII[) . Eq. (TT21 implies 



GK eq ( i ,l,w) = 2iG n (i-j,0)J 3 - 



and 



E8{uj)G A {k- 1,0) 
(14) 



— 04 T 



E5(u) 



(15) 



For small E, hs is just two powers of E times an inte- 
gral representation of the derivative of the delta function 
so terms involving hs would seem to give only contri- 
butions of higher order in E. But as shown in Ref. [ID 
and discussed in more detail below, the singularity in the 
T=0 equilibrium distribution function means that the hs 
terms give a contribution of order \E\ to the noise. This 
issue will be further discussed in Section IIVI 



D. Path Integral Representation 

We write the model as a path integral on the Keldysh 
two-time contour^^! For conducting systems which 
have magnetic ground states we expect that part of 
the interaction term may be decoupled by a Hubbard- 
Stratonovich transformation™ while the remainder can 
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be treated by Fermi liquid theory. The important part of 
the Hubbard-Stratonovich field is a vector in spin space 
rh(r,t) ■ a which represents the local magnetization. 

After integrating out the fermions we may express the 
model as a path integral 



[{m(r,t)}]+iS int [{m(r,t)}] 



T>me 



(16) 



In Eq. (|16p we consider time as running along the 
full Keldysh contour (including forward and backward- 
moving parts and with 'Keldysh' time ordering assumed). 
The term Si n t[{m(r, t)}} includes other terms in the ener- 
getics of the local spin polarization, including the terms 
arising from the Hubbard-Stratonovich transformation 
and the magnetic dipole energy. Sp , the logarithm of 
the fermion determinant, is 

S F = -iTrln [G^ 1 - mn,(f) • oV 3 - T 3 S to tT 3 ] (17) 

Sf depends on Gq , which is a matrix in space, time 
and Keldysh contour indices and expresses the effects of 
band structure and of Fermi liquid renormalizations. The 
matrix f acts in Keldysh space. 

Following Schula2£ we now rotate the spin quantization 
axis at each space-Keldysh time point so that it is parallel 
to n. The set of space-time rotation matrices which do 
this arc R„ defined so that Rt n <rRfi = a z . Writing the 
spin orientation vector n(r, t) = cosOz + sinOu with u a 
vector in the x — y plane, it is convenient to represent R 
in terms of defined by 



O = I cos—z 
2 



r = n,i 



stn—u 
2 



x ■ a 



(18) 



(19) 



We obtain 



iTr In 



Rt G^R fl - m<i z - Rt T 3 £ tot T 3 R A 



(20) 

The spin rotations correspond to an SU(2) gauge 
transformation of the theory. If all of the terms in the 
underlying Hamiltonian are invariant under rotations of 
spin in the system of interest then the trace must be 
independent of the choice of spin quantization axis, cor- 
responding to a gauge invariance of the model. Terms 
such as the local magnetization ran • a or sources (not 
explicitly written here) will break the SU(2) invariance 
and lead to dependence on the spin directions. The reser- 
voir which provides spin relaxation also breaks the SU (2) 
gauge invariance because as we have defined it the rota- 
tion does not act on these degrees of freedom. Thus for 
example a time dependent spin rotation corresponds to 
a spin-dependent phase difference between system and 
reservoir, which could drive spin currents between the 
system of interest and the reservoir. A similar situation 



obtains in the theory of dissipative effects in supercon- 
ductors^^ 

From Eq. (|2"0")) , using RtR^ = 1, and defining 
the rotated frame Greens function G via G 1 = 
(G^ 1 — ma z — t 3 Et 3 ) we obtain 



Here 



and 



S F = -iTr In [G" 1 + SH - <5£] 
<5H = Rt G 1 R„ — G 1 
<5£ = Rtr 3 Sr 3 R ft - r 3 Sr 3 



(21) 
(22) 
(23) 



Assuming the magnetization varies slowly in space and 
expanding to second order in rotations gives 



with 



S — So + Si + S2 

50 = -iTrln [G _1 ] + S int 

51 = -iTr [G (<5H — <5£)] 

5 2 = ^Tr [G (SH. - SY,) G (5H - S'E)} 



(24) 

(25) 
(26) 

(27) 
-S 

and retain terms up to first order in departures from 
equilibrium. Considering first the equilibrium terms we 
observe that our assumptions of non-magnetic scatter- 
ing, slow variation of the magnetization and relaxation 
of quasiparticle distribution to the local magnetization 
imply that only the terms Ejv-c arising from spin non- 
conserving quantities contribute to 5H eq : 



At this point we write G, £ — > G 



eq 



Rt T 3 £^_ c T 3 R fi - t 3 5^_ c t 3 (28) 



Thus for the equilibrium terms we have 



8S eq = Sl q + S 



with 



Sl q = -iTr[G eq (SH-S-E e ^ c )] 



? e 9 _ 



-Tr 



G eq (SH - ST, eq _ c ) 



:G eq (SH-S-E eq _ c 



(29) 
(30) 

(31) 



Here S^ 9 will be seen to give the basic Berry phase mag- 
netization dynamics while S^ gives the spin stiffness. 

Turning now to the nonequilibrium terms, we work to 
leading order in both departures from equilibrium and ro- 
tations, neglecting the unimportant change in spin stiff- 
ness due to current flow. We obtain 

s ne g = _ lTr [ Gne9 ^ H _ S ^_ c )] + iTr- [G eq 5V neq ] 

(32) 
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The rest of this paper is based on an evaluation of 
Eqs. (|5D|). (|3Tj). (|3"2j) in the presence of spin tex- 

tures and externally imposed currents and in the limit of 
slow variation of spin direction. We perform an expan- 
sion in the amplitude of the quantum components of the 
rotation matrices. The leading terms in the semiclassical 
expansion give rise to the classical equations of motion 
for spin fluctuations and topological defects in the pres- 
ence of an applied current. The next-to-leading terms 
lead to a current-induced effective temperature. 



III. DRIVE AND DYNAMICS 

A. Overview 

In this section we present the leading terms in a semi- 
classical expansion. These terms give rise to the classical 
equations of motion for spin fluctuations and topologi- 
cal defects in the presence of an applied current. They 
are most easily obtained by going to the Larkin basis in 
contour spaced introducing classical (C) and quantum 
(Q) combinations of the time-ordered R2 and anti-time 
ordered Ri Keldysh contour fields via 



where the trace is over space, time and spin degrees of 
freedom. 



Now 



Go'~Hn - r 2 ;h - t 2 ) = 6(h - t 2 ) (id t - 



and G A ' 1 (ri 



n v 1 ~h) = C-id t - tij^J S(ti-t 2 ) are 

short ranged in space and time so we may perform a 
gradient expansion of Eq. (12SD. Because R Q (1)R C (1) + 

R^(1)R Q (1) = (R|(1)R 2 (1) - R t 1 (l)Ri(l))/2 vanishes, 
the leading terms are the first derivatives. The space 
derivative term vanishes because in equilibrium G K is an 
even function of its space arguments. Thus we obtain 



Rc 
Rq 



2(R 2 + Ri) 
^(R 2 -Ri) 



(33) 
(34) 



and expanding to leading order in Rq. We evaluate the 
terms in order. 



eq,G0 



-iTr 



G* ({-id t Rl}Il Q + R Q {id t R c } N 

(36) 

-iG K (i, i; t, t) = p + ma z with p a scalar re- 
lated to the charge density and m the magnitude of the 
magnetization on site i at time t. z is the (arbitrarily 
chosen) spin quantization axis. 



We first consider the terms involving SH . Expanding 
and taking the trace over contour indices gives 



ST 



eq,G _ 



-iTr 



G K {\- 2)Rt 7 (2)G^" 1 (2 -l)Rg(l) 



G^(l-2)R t o (2)G «'- 1 (2-l)R c (l) 



(35) 



We next consider the term in «S^ 9 arising from the equi- 
librium self energy. Expressing the result in the Larkin 
basis, taking the trace and retaining only the terms with 
one quantum field gives (here we leave the trace over spin, 
space and time labels implicit, denote spin matrix struc- 
ture in bold and collapse the space and time indices): 



Sl eq = iTr G fl (l-2)R^(2)E^ C (2-1)R Q (1) + G K (1-2)R^(2)I]^ C (2-1)R Q (1) 

+ G A (l - 2)R Q (2)££_ C (2 - 1)R C (1) + G*(l - 2)R Q (2)£#_ C (2 - 1)R C (1) 

I 



(37) 



Our assumption of spin-independent reservoirs means 
that E R ' A ' K are scalars in spin space so may be factored 
out, leaving 



Sf >eg = i D CQ (l - 2)R t c (2)R Q (l) + (QhC) (38) 
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Here, 

B CQ = G«(l - 2)£*_ c (2 - 1) + G K (1 - 2)S^_ C (2 - 1) 

(39) 

and 

D QC = G A (l - 2)£*_ c (2 - 1) + G A '(1 - 2)S«_ C (2 - 1) 

(40) 

Since G R I A and S^^ c are short ranged in space and 
time we can make a gradient expansion of the prod- 
uct of rotation matrices, writing R^(2) = R^(l) + 
(2 - l)<9R f c (l) + 0(1 - 2) 2 and R c (l) = Rc(2) + (1 - 
2)dR^(2) + 0(1 - 2) 2 . We then change labels (1 <-> 2) 
in the second term and obtain 



d2 



D CQ (l-2)R t c (l)R Q (l) 



+D QC (1-2)R Q (1)R C (1) 



(41) 



Since f d2D C Q(l-2) = J d2D QC (l-2), the leading term 
vanishes while the even parity of the equilibrium state 
means the first space derivative vanishes. So the leading 
term is the first time derivative and incorporating the i 
into the coefficients we obtain 



with 



Tr 



D 



CQ 



D 



QC 



Dc Q a t R t c (l)R Q (l) + Dq C R^1)5*R c (1 

(42) 



/ / d2(t 2 -t 1 )T>c Q (f 1 -f 2 ,t 1 -t 2 ) (43) 
d2{t 1 -t 2 )T> QC (f 1 -r 2 ,t 1 -t 2 ) (44) 



The coefficients are most conveniently evaluated in 
Fourier space as D C g = c?Dcq(0) /90|n=o and Dqc = 
-9Dgc(^)/9fi|n = o- Introducing 



G 
A 



_(G^(l-2) + G«(l-2)) (45) 

I(G A (l-2)-G*(l-2)) (46) 

~(s£_ c (l-2) + £*_ c (l-2)) (47) 

I(E£_ c (l-2)-££_ c (l-2)) (48) 



and noting that in equilibrium we have G K (k,u>) — 
—2iA(k,Lu)h(uj) and £^_ c (fc,u;) = — 2iTjsr-c(k, u)h{w) 
with /i the equilibrium distribution function we find 



(49) 
(50) 



D CQ {tt) = Di(Sl)-iD 2 (Sl) 
D QC (n) = -D 1 (Sl)-iD 2 (Sl) 

with the real functions Di and D 2 given by 

Di(O) = 2 J{dkduj)A(k,uj + n)r N -c{k 7 uj) 

x (h(uj + tt) - h(uj)) (51) 



and 

D 2 (fi) = 2 J{dkduj) G'(k,u} + ir)T N - C (k,Lj)h{u>) 

+A(k,u + n)Z N _ c (k,u)h(uj + n) (52) 

Thus, L> C q(^ = 0) = L>qc(^ = 0) as noted above. We 
therefore obtain 

S l^ q = T,j dt 

i 

Dberry (-M^R^i, t)K Q (i, t) + Rg(l, t)ld t Rc(«, *)) 

+D dlss (a t Rj,(i,t)R Q (i,i) +Rj } (t,t)ftRo(*,*))(53) 



with 



D 



diss — 



7m 



^ J (dfc)A(fc, cj = 0)rAr-c(fc, w = o) 



(54) 



and, using the Kramers-Kronig relation and shifting uj — > 
w — f2 appropriately, 

Dfc errj , — ~2 y (duidx)(dk) 
A(k,x)T N -c(k,oj) (h(u) - h(x)) 

/ \2 ( 55 

{UJ — x) z 

The T)b erry term appears as a renormalization of Eq. ([36]) 
while the D<jj SS term parametrizes effects of dissipation. 
In the rotated basis both terms are diagonal matrices 
in spin space and are explicitly proportional to the spin 
relaxation rate Tn^c- Additional terms would appear if 
the reservoir was magnetic. 

We have written Eq. (f5~4"f for the general case. Our spe- 
cific results are presented for a momentum-independent 
spin relaxation rate. For this case the Tn-c may be re- 
moved from the integral and we identify f(dk)A(k, u> = 
0) as the spin-dependent Fermi level density of states 
N a . In a simple model with k 2 /2m dispersion, N a = 
mkp^/ (2ir 2 ). Thus we estimate that the diagonal ma- 
trix elements D" are 



_3n a T N -c 

diss 



7T E Fm 



(56) 



with Ep a = k Fa /2m the Fermi energy for spin direction 
a. 



C. S° q 



Expanding , Eq. (|3"Tj) . gives the standard spin- 
stiffness terms, renormalized weakly by the coupling to 
the reservoir. These may be written (reverting here to 
the physical spin coordinates) 

S stiffness = ips J dxdt h q ■ V 2 n cZ ^(x,i) (57) 

where n c ; t ±(x,t) is the classical component of n perpen- 
dicular to the local magnetization direction and ps is the 
spin stiffness, proportional to the local magnetization. 
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D. Nonequilibrium terms 



We first consider the first term in Eq. (|32[) . Expanding 
as in the previous section we find terms of the form of 
Eq. (I35[) but with G^ replaced by the nonequilibrium 
Keldysh function and terms of the form of Eq. (|37|) with 
the G R ' A terms removed and the G K replaced by the 
nonequilibrium Keldysh function. The terms of the form 
of Eq. (ESI) are 



s neq,G = _ iTr 



G£,(l - 2)R^(2)G^- i (2 - l)Rc(l) 



G* (1 - 2)Rt 7 (2)G^" 1 (2 - 1)R Q (1) 



R-l, 



Since the nonequilibrium G^ eQ is an odd function of 



its space coordinates but an even function of time the 
only non- vanishing contribution comes from the hopping 
terms in SR. We obtain 



S. 



neq,Go 



iTr 



G ne „(*, ji fyj—i 



x (Rk(j,i)R Q (M)+Rk(j } i)R c (i,*) 



(59) 



Writing R^(j) = (fj - r l ) a V a Rj 7 (i) and R c {i) = (n - 
fj-)°V°Rc(i)) changing i to j in the second term and 
using Eqs. (JHJ), ([7J we obtain 



-iTr 



r (w a nl(i,t)R Q (i,t) 



R' (z,t)V a R c (z,i) 



(60) 



We next consider the terms arising from S eg . These 
are of the form of Eq. (T3T|) but with the functions 
D QC/CQ replaced by 



and 



Ecq = G£,(1-2)E£_ C (2-1) (61) 



^QC = G*,(1-2)E*_ C (2-1) (62) 



Making the gradient expansion as above gives 



Tr 
with 



E 



£ Q V a R^(l)R Q (l) + E^ C R^(1)V Q R C (1) (63) 



CQ -= i j d2(r 2 - nY^cQin - r 2 ,h - t 2 ) (64) 



E QC = i J d2(fi -r 2 ) a E QC (fi -f 2 ,ii -t 2 ) (65) 

In the limit of momentum-independent (space local) spin 
relaxation rates which we consider here these terms van- 
ish, but in general there will be a contribution propor- 
tional to dk'F'N-C- 



In the second term in Eq. (|32[) . use of the steps lead- 
ing to Eq. (|37l) along with the fact that T, neq has only a 
Keldysh component and is proportional to the unit ten- 
sor in spin space gives (leaving the sum over coordinates 
implicit) 

Sf, neq = i F CQ (l-2)R t c (2)R Q (l) 

+ zF QC (l-2)R t Q (2)R c (l) (66) 



with 



(58) and 



F CQ = G R (1 - 2)E£,(2 - 1) (67) 



F QC — G A (l-2)E r f e „(2-l) (68) 



Again making the gradient expansion and incorporating 
the i into the coefficients we find 

SSl neq ^ = F£ Q V Q R^(1)R Q (1) + F% c R t Q (l)V Q R c (l) 

(69) 

The explicit form of E^, (Eq. (fT5|) ) gives, in a mixed 
position- frequency representation 

f cq = \ E ~ * u = [fa ■ 3) ( 7 °) 

3 

F QC = ~ E G ^ " 3, " = Ojjti-i^ (f* • (71) 

Here r*y is the vector connecting site i to site j and in the 
QC term we have r*- not r* 4 because we put the gradient 
on the Re- In Fourier space we have 



CQ 



(dk)G R (k,Lo = 0) 



d 2 e k 
dk a dk b 

,2 



E b (72) 



FJ C = 4/ W G^ = 0,|S- E » <™> 
Our final expression therefore is 

Sf,neq = E/ / ^ 

iF} (-V a R^,(i, t)Ro(i, t) + Rq(«, t)V a R c (i, i)) 
+F^ SS (v°R^, t)Ro(i, t) + Rt,(i, t)VRc(i, *)' 



(74) 

with the F diagonal matrices in spin space with up and 
down spin components (here we use the resistivity to ex- 
press the electric field in terms of the current) 

F?_ = -- / (dk)TZeG R (k,uj = )-fg-p CT 475) 



J,<T 



dk a dk b 



(76) 
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It will be useful for our further discussion to estimate 
the magnitudes of F a in a simple model with k 2 /2m dis- 
persion. In this case d 2 £k/ dk a dkb = S ab /m and the esti- 
mates that led to Eq. (|5r?|) along with Eq. (JS| give 



diss, a 



3ri 



2-kE f 



-J a 



(77) 



S c uss =J2 Tr [ D m{ r q (^Rc) + (fyRk) Rq} 



(84) 

S stiffness is given by Eq. (J5TJ) and S noise by Eq. (79]). 
Here /i=0 denotes the time direction and fi=a a space 
direction. 
We have 



IV. SECOND ORDER IN QUANTUM FIELDS: 
NOISE TERMS 

Expanding the action to second order in quantum fields 
yields terms which stabilize the integral over quantum 
fields and, among other effects, introduce an effective 
temperature into the problem. A general analysis of 
these terms is somewhat involved; we focus here on ob- 
taining the qualitative behavior at T—0 and to leading 
order in the gradient expansion. In this case the equilib- 
rium terms do not contribute and expanding Eq. (132[) to 
second order in Rq yields 



iTr 



(78) 



For simplicity we assume that the spin-nonconserving 
scattering rate is small compared to the total scatter- 
ing rate, so that the first term is dominant. We then 
obtain 

Snoise =*T*J2 I dtTr [NR^, t)T tot R Q (i, t)l (79) 



N is a diagonal matrix with entries N a = 
^ J(dk)\cos(j)k\A a (k,uj=0), <f>k is the angle between 

Vk and E, T is a diagonal matrix with entries yT^*T*°* 
with 



. . r*°* + r* ot 

ptot _J 4- 

2 

and the effective temperature is 

\E\\v F \ 



(80) 



(81) 



As in Ref. [l3|, the effective temperature is given by the 
voltage drop over the distance traveled between inelastic 
scattering events. 



(85) 
(86) 
(87) 



Bo — m<J z + T)berry 

B a = 3 a + F a j 

E>0 = 'Ddiss 

r\ T?a 

'-'a — x diss 

Bold-face coefficients denote matrices in spin space; the 
spin structure arises from spin dependence of the density 
of states and scattering rate in the presence of a non- 
vanishing local magnetization. 

We now use Eq. (TT~5T) to express Eq. ([8"3]l in terms of the 
classical and quantum components of the unit vectors f2. 
We have 



Rn = flnl + iz x tin ■ a 



and 



d u R, 



<9R, 



The action is thus of the form 



ibr^b 



-iT*H,QN a fig + S s tiffness 



with 



jab 



iTr 



B, 



D, 



'dn c dR c 



dK c dR c 



on a an b on b an a 



c gft-c 



an a dn b dn b an a 



and 



N xx = N yy 



iV T r*°* + N\ T tot 



N t T\ ot + NiT\ ot 



(89) 
(90) 

(91) 
(92) 

(93) 

(94) 
(95) 



V. SUMMARY OF ACTION 

Collecting all of the terms we obtain an action in terms 
of the classical and quantum rotation fields of the form 

S F [{Rt,R}] =S B + S dlss + S 

stiffness "T" S no i se (82) 

with 



B„{Rk (id n K c ) ~ (id^c) Rq} 



(83) 



Evaluating the derivatives explicitly and noting that 
J ab and D ab have terms ~ 1 and ~ a z we find 



rxy _ _ ryx _ r, rr 
• J a ° u. ~ a 



with all other elements of while 



Df = 



ab 



(96) 
(97) 



We thus have 

S B + S dlss = ti Q ■ (-z x (J^) ti c + {V^) tic) (98) 
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with [i = 0, a labeling space-time components and 



Jo 



2Tr \a z T> 



berry J 



J a = 2Tr[a z {J a +F 
V Q = 2Tr[T) dlss ] 
V a = 2Tr[F% ss ] 

The estimates provided above imply 



Jo 
J 

V 
V 



Jm 



mj 



E F 



J 



(99) 
(100) 
(101) 
(102) 

(103) 
(104) 

(105) 
(106) 



The coefficient of VLq in Eq. (|9"5|) gives the classi- 
cal equation of motion for an isolated spin in the pres- 
ence of dissipation and current drive. It is consistent 
with previous derivations based on the phenomenological 
Landau-Lifshitz-Gilbert equation; 16 ' 17 ! 22 and with micro- 
scopic derivations of the Landau-Lifshitz equations in the 
presence of current i 23 ' 25 Here the terms proportional to 
J are non-dissipative terms giving rise to the adiabatic 
spin torque and arising from the usual Berry phase, as 
modified by an applied current ) 18 ' 21 ' 29 while the terms 
involving T> arise from dissipative processes and give rise 
to the nonadiabatic spin torque^ 2 - 

The Galilean transformation 



with 



r — Uvt 



J_ 

Jo 



(107) 



(108) 



shifts to a frame co-moving with the spin current. In 
this reference frame the spatial derivative terms in Sb 
vanish. The coefficient of the spatial derivative in the 
combination of Sdiss becomes 



"Do (Un 



with normal velocity 



V 
V~o 



(109) 



(110) 



In a Galilean invariant theory this combination vanishes. 
The easiest way to achieve Galilean invariance is to re- 
move the scattering, so 2^=0, however one may also 
imagine a situation in which u n — u s . However, the esti- 
mates given above show that in the general case, in the 
co-moving frame the coefficient of the spatial derivative 
in the T> term is of order of the difference between the 
transport and spin relaxation rates (measured relative to 
the Fermi energy) 



r N - c - r tot 

E F 



J 



(111) 



Because we generically expect a spin relaxation rate much 
less than the transport relaxation rate, Galilean invari- 
ance is strongly broken in general. 



VI. EQUATION OF MOTION FOR A 
TOPOLOGICAL DEFECT 

A. Overview 

To obtain the equation of motion of a vortex we imag- 
ine that spin wave excitations can be integrated out so 
that the spin orientation n is determined by topological 
defects at positions A J '(t). We further suppose that the 
topological defects move slowly enough that the spins are 
relaxed to the values appropriate to the instantaneous po- 
sition of the defects. We have (repeated indices summed) 



n a = 



(112) 



Then we may rewrite Eq. 



as 



ia ijab 



+iJ2 x Q,M sahx Q,j( t ) ( 113 ) 



Here the noise term 



r?ab 



5Q? 



SXg./ ■ 5X^ 



(114) 



The generalized force T and Liouville operator C are for- 
mally given by 



5Q 



5X. 



— ■ (v ■ vn) - V 



dU(X Ctl - X Clj ) 



(115) 



£#(t) = / (dr)Joz 



sn 



Do 



sxk; 



(116) 

Here f2, z denote unit vectors in spin space. They 
depend on all of the vortex positions, although this de- 
pendence is not explicitly notated. The V are vectors in 
position space. Except for the terms in parentheses in 
Eq. 



(jl 15|) the dot and cross product structure refers to 
the spin directions. J is the spin current produced by an 
applied electric field. The direction of current flow (indi- 
cated by the vector symbol) is along the applied electric 
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field while the spin direction is parallel to the local mag- 
netization. V is also directed along the applied electric 
field. U is the interdefect potential arising from the spin 
stiffness and dipolar forces. 

The equation of motion for vortex i is obtained in the 
usual way2£ by using a Lagrange multiplier to integrate 
out the quantum fields, leading to 

rjCJoZ xX + I D V X t = v CzxJ + I D V + ^ F tJ + Q 

with Fij — —VijU(Xi — Xj), n = ±1 the sign of the 
vorticity and Q a fluctuating field which on scales longer 
than the mean free path may be taken to be delta corre- 
lated, with 

(G(t)$V)=2B$6(t-t) (118) 

Eq. (|117[) is consistent with derivations based on the phe- 
nomcnological Landau-Lifshitz-Gilbert equation ] 16 i 17 

In the remainder of this section we examine separately 
the drive terms (those arising due to the presence of an 
imposed current and visible in the equation of motion for 
a single topological defect) and the intervortex interac- 
tion terms. 



B. Drive terms 

We examine Eq. (|116[) for a two dimensional system 
containing a single vortex, at a position X(t) = X(t)e x + 
Y(t)e y . In this case Q -> Q(f - X{t)) so SQ/SX a = 
— V°fl and Eqs. (|115[) . (|116|) involve two quantities: 

h = J d 2 rM{r)(v x n x Vyh y -v x h y v y n x yii%) 

Id = J d 2 r\7 a {l ■ V 6 (120) 

with Ai proportional to the magnetization a function to 
be discussed below. Ij and Id depend on the structure of 
the vortex. We consider here an easy-plane system with 
spins lying preferentially in the x — y plane and choose 
the origin of coordinates to be the position of the vortex. 



Far from the vortex, M. is a constant proportional to 
the local magnetization and the spin orientation vector is 
h x — cosi](f>, hy = sinr](f) and h z = 0, implying £l z = 

and (fijc, £l y ) = -^={cosr)(j), sinr)(f>). Within a distance £ of 
the vortex the structure changes. The integral defining 
Id is dominated by the far region. Substituting we have 

I a D h = m f8 ab lnj + ... = I D 8 ab (121) 

with the ellipsis denoting terms of order 1. In a system 
with a single vortex the long distance cutoff is the sys- 
tem size; in a system with equal densities of vortices and 
antivortices the cutoff is the mean intervortex spacing. 
This logarithm has been previously notedi^ 

The noise term may be evaluated similarly. The lead- 
ing behavior is also logarithmic and arises from the case 
where i and j refer to the same vortex; we estimate 



S a = I D N aa T* (122) 

and S z of order unity. Here T* is given in Eq. (|81[) . 
Note that in a fully polarized (half-metallic) magnet, 
N xx = and the dissipation vanishes in this approxi- 
mation (higher order processes will lead to a non-zero, 
but much smaller, noise). 

Turning now to Ij we see that the contributions from 
the far region vanish so Ij is determined by the properties 
of the vortex core, which depends upon details. Two 
simplifying limits are possible. In the extreme soft spin, 
high anisotropy limit, the core is defined by a vanishing 
of the local magnetization amplitude (represented here 
by M. — > 0) while the spin vector remains in the easy 
plane. In this case lj=0. On the other hand, in the hard 
spin, weak anisotropy limit the spin amplitude remains 
constant, and the spin direction rotates out of the plane, 
becoming parallel (or anti-parallel) to z. In a general 
situation both effects occur, with the relative importance 
controlled by the ratio of the magnetic anisotropy energy 
to the amplitude energy. 

Going to polar (r — <\>) coordinates we have Cl x = 
sin^r-cosrjcj) and fl y = sin^r^- sinqcf) so 



j2 » j, n f / dB 6 2 , rj . e . « ,,\ / dB 9 . , , t? . 9 2 ,, 

a rA4(r)< cos— cos rim H — sm— sin rim) cos— sin nm -\ — sin— cos rim) 

y '\ \2dr 2 ' r 2 1 ' ) \2dr 2 ' r 2 ' ; 

/ d6 9 , . , n 9 \ ( dQ 6 f] . 9 , A 1 , . 
——cos—cosn(psinrj(f) sin—cos(psinrj(f)) I ——cos—cosrj<psinn(p sin—cosrjcpsinrjcp) I > (123) 



The terms which survive are At r = oo cosO — and at r = cosO = ±1 according 

r fja %n f dcosO as ^ e s P m rotates parallel or antiparallel to the z axis. 

I j =rj d 2 rM(r) — —sinO = — / drM(r) 



Ardr 2 / dr 

(124) 
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Thus if M. is constant Ij is a topological invariant, pro- 
portional to the vorticity and to the sign of the rotation, 
but variation of the magnetization amplitude near the 
core means the amplitude is non-universal. In the ex- 
treme hard axis limit the spin remains in-plane over the 
whole region where M 7^ and Ij vanishes. In general 
we can write 



Ij = vC 



(125) 



with C of order 1 carrying the sign of the spin orientation 
in the vortex core and the magnitude factors associated 
with the variation of M. near the core. 
Thus we finally obtain for the action 



5' 



: I ' dtY,T*N xx i D x Ql - X Qi 

J i 

f J2i D x Qi - (v-v q x 



-X ( >, ■ { n C: x \ J-J Q k Ci . 

du{x c ,i - x ( , ) 

j 



E^-E 



dX c , 



With the help of Eq. (fT08| we identify 



F 



M 



= vCJaz x (u s - X^j 



(126) 



(127) 



as the Magnus force acting on a vortex. We see that it 
has a sign determined by the combination of the vortic- 
ity and the constant C which describes the out of plane 
direction of spin rotation in the vortex core. Compar- 
ison to the usual theory of neutral superfluids (see e.g. 
Eq. (2.1) of Ref. Hil) allows us to identify J with the 
combination 2irps /m (essentially 2ix times the density of 
superfluid particles). On the other hand, to make con- 
tact with results obtained for superconductors it is more 
useful to identify J = J§u s with the supercurrent and to 

note that the combination JqXq gives rise to the vortex 
Hall angle, which in superconductors is typically of order 
T c /Ef and is in most cases negligible. 

Proceeding similarly and using Eq. (|110l) we see that 



F D = I D V Q (u n -X 



(128) 



represents the physics of dissipation. A very similar term 
occurs in superconductors and in the theory of neutral 
superfluid films (see e.g. Eq. (2.2) of Ref.Ql)," although in 
this case physics not relevant here produces an additional 

term proportional to z x (u n — X). 

In neutral superfluids and superconductors, u n — X 
represents the velocity of the vortex relative to the source 
of dissipation (the substrate, in the neutral superfluid 
case, or the normal-fluid excitations in a two-fluid model 
of the superconducting case). The relative magnitudes 



of u s and u n is an important issue in the theory of su- 
perfluids shaken by motion of their container, but except 
perhaps in the ultra-clean limit u n may be neglected in 
superconductors. However, in the magnetic case our es- 
timates indicate that u s /u n ~ r 'n-c /T tot is typically 
much less than unity. 

The relative magnitudes of the dissipative (IdF>q) and 
non-dissipative (j7o) terms sets the vortex Hall angle. In 
typical superconducting films the dissipative terms are 
dominant, but in the magnetic situation of relevance here 
the order of magnitude of the dissipative coefficient is set 
by the spin relaxation rate Tn_c which is likely to be 
small relative to the other terms in the vortex equation 
of motion. 

For a single vortex (so F=0) and in the absence of 
noise we may write the solution of Eq. (|117l) explicitly as 



dX 
~dt 



r)CJ I D V 



(i D v y + (riCJ y 



z x (u s - u n ) 



(IdVoY 



(i D v y + ( v cj Q ) 



(u s -u n ) + u s (129) 



We observe that if either the spin relaxation rate (and 
hence T>q) vanishes or if u s = u„, Eq. (|129l) implies, 
X = u s , so that the vortex simply drifts with the ap- 
plied current. Thus in order to get motion transverse to 
current flow one needs both non-zero spin relaxation and 
a breaking of Galilean invariance u s 7^ u n . In our model 
both of these effects arise from the presence of the spin 
reservoir. 



Eq. (|129[) is very similar to the equation of motion for a 
vortex in a neutral superfluid given in Eq. (2.3) in RefQl 
which we rewrite below, 



dR 
~dt 



B 



Hsl y{(2nh/m)p° s - B') 2 + B 2 
+C (v n -v s )+ v : 



I x (v n - v s ) 
(130) 



Here v s is the superfluid velocity, while v n is the veloc- 
ity of an underlying substrate which is also responsible 
for the dissipation and n = ±1 denotes the vorticity. B 
is a drag coefficient arising from coupling of the super- 
fluid to the substrate, and B' is a drag coefficient aris- 
ing from coupling of ripplon excitations to the curl of 
the superfluid velocity which does not appear in the case 
considered here. When B' = 0, the coefficient C defined 
in Ref. 19 equals B 2 /((2Trh/m)p t) s ) 2 + B 2 ). Identifying 
B with IdT^o and {2irh/m)p° s with CJq we see that the 
expressions become identical. 

The important term for our subsequent considerations 
is the first one on the rhs of Eq. f|129[) , which shows that 
vortices have a component of motion which is transverse 
to an applied current with sign determined by the combi- 
nation -qC. The estimates presented in the previous sec- 
tion show that 2?o is of the order of the ratio of the spin 
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relaxation rate Tn-c to the fermi energy Ep while T> is 
of the order of the ratio of the total scattering rate T tot 
to the fermi energy. Jq is proportional to the magnetiza- 
tion while the spin current J is proportional to m times 
the total current J. Thus the transverse component Xp 
of vortex motion is roughly 



- (nCz x J) ■ 



ml D - 



Id- 



E f 



- (r/mC) 2 



(131) 



Except very near to a magnetic transition or in the 
case of very strong spin relaxation the first term in the 
denominator is negligible and we have 



X, 



z x J 



\mri 2 C 2 



Interdefect Interaction 



(132) 



The force Tij between two vortices at positions Xi and 



Xj is given by 



(133) 



It has contributions from the dipole force and from the 
usual spin stiffness terms. Writing U = Uss + Udi P , the 
interaction arising from the spin stiffness is 



Uss 



In 



Xl2 



(134) 



with ps the spin stiffness per unit area. 

In magnetic systems, unlike in superconductors, the 
dipole interaction may also be important. For a general 
spin configuration the dipolar energy may be written as 

U dip = (gn B ) 2 / d 3 xd 3 x' (135) 

J \X X \ 

with g the electron g- factor. Specializing to a two dimen- 
sional film of thickness d and considering the case of a 
vortex and an antivortex separated by the vector X12 we 
obtain^ 



Udip = Vi^EdipiX^ ■ m) 



(136) 



Here the dependence of X follows by power counting and 
Edip = {gpB) 2 d 3 Cdi P oie(Xi2-m) is the dipolar energy as- 
sociated with magnetization in a volume set by the film 
thickness d, multiplied by a coefficient which depends 
on the relative orientation of the vortex/antivortex sep- 
aration and the magnetization. We have evaluated the 
constant C 'dipole numerically, finding C 'dipole ~ HO for 
a vortex/antivortex pair separated along the direction 
defined by the magnetization and Cdipoie ~ 35 for a vor- 
tex/antivortex pair separated perpendicular to the direc- 
tion defined by the magnetization. 



VII. CURRENT-DRIVEN VORTEX 
UNBINDING TRANSITION 

The previous section showed that vortices of opposite 
vorticity tend to be driven apart by an applied spin cur- 
rent but are pulled together by the intervortex force. The 
competition between these two effects leads to a current 
driven vortex unbinding transition. To investigate this 
possibility in more detail we consider the equation of mo- 
tion of a single vortex-antivortex pair, obtained from the 
solutions of Eq. (|117[) for two vortices of opposite vor- 
ticity (more precisely opposite sign of r]C, and in what 
follows we take \r/\ = 1). We find that the equation for 
intervortex distance X±2 can be written 



with 



X ia =A[F-£)+C 



A = 



2I D V 



(I D V ) 2 + (CJo) 2 



(137) 



(138) 



The force term coupling the two vortices is given by 
Eq. (|133p . The current-induced force acting to separate 
the two vortices is 



(139) 



We observe again that the net force vanishes if the 
Galilean-invariance condition JqD = T>$J is satisfied. 
£ is a noise field with correlations now given by 



<C«(0 = 2AT eff S lJ S ab S(t - (140) 



with 



T eff = 



Vn 



(141) 



Eqs. (|137|) . (|136|) show that if the applied current is 
large enough that 



V a 



(142) 



then at sufficiently long distances it becomes energeti- 
cally favorable to separate vortex/antivortex pairs. 

To calculate the probability of formation of vor- 
tex/antivortex pairs we follow Ref. [l^ and recast 
Eq. (|137[) as a Fokker-Planck equation for the probabil- 
ity P per unit area of finding a vortex-antivortex pair at 
separation r: 

d t P(r, t) = -V • J P (143) 
Here, the probability current Jp is given by 

J P = -AT eff e~^VP{r,t)e^ (144) 
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and effective potential 



U 



E, 



dip " 



np s \r\ 
2 d 



Ex ■ 



(145) 



In the limit of relatively weak applied currents, 
Eq. (|143p may be analyzed by saddle-point method i 19 ' 31 
We take the electrical current to be applied in the y direc- 
tion, so the vortex/antivortex pair separates along x and 
choose x = y — as the location of the saddle point in 
U with x — ¥ oo corresponding to large vortex/antivortex 
separation. We assume that the probability current den- 
sity is only large near the saddle point of U and seek a 
steady-state solution with J v p — 0. The pair creation rate 
per unit area R is then given in terms of the y integral 
of the x component of the probability current a s 19 ' 31 



R = J dyJp{x = Q,y) 



(146) 



To determine R we observe that the condition that 
is fulfilled if P has a y dependence which com- 
pensates the y dependence of U, i.e. 



j p — 



(147) 



For x — s- oo the density of vortex/antivortex pairs is 
very low, so / — > 0, whereas for intervortex separations 
much less than the saddle point value we expect the prob- 
ability distribution to be the pseudoequilibrium one cor- 
responding to temperature T e // so f(x « 0) — > 1. The 
steady state condition is fulfilled if Jf, depends only on 
V, i-e. 



T x 



$(y) = -AT eff e T -ff d x f(x) 



(148) 



Multiplying both sides of Eq. (|148[) by e T ?ff , inte- 
grating over x using the boundary conditions on / and 
rearranging then gives 



AT, 



eff 



U(* ,y) 

f°° dx'e^TT 



(149) 



with xi ow sufficient far from the saddle point that 

To implement these formulae we now make the saddle 
point approximation explicit. The saddle point location 
r c is 



Near the saddle point we write 



U(Sx,Sy) _ 5y 2 8x 
T eff = U+ 2 y *~2x 2 



with 



nps r c 
-in- 



2T eff d 2T eff 

-2 _ Edip nps 

T eff dr c 2T eff r 2 
-2 _ ^ps 



IT. 



eff 



(150) 

(151) 

(152) 

(153) 
(154) 



Thus Eq. (I148P becomes 



J x P (y) 



AT, 



(155) 



R = 



1 f2e8Sd\ V ° 2 "* iT * 



x c v2tt 

so, performing the y integral and restoring units 
2I D N XX T 

dse 

(156) 

Because the vortex annihilation rate scales as the 
square of the vortex density and must balance the cre- 
ation rate, we obtain 



(I D V ) 2 + (CJo) 2 Ji,Kur V vps ) 



Xly 



[288d\ D °zw**r* 



(157) 



VIII. CONCLUSIONS 

In this paper we have shown that an applied cur- 
rent can drive a nonequilibrium topological (defect- 
unbinding) transition in a continuous-symmetry magnet. 
We derived the appropriate action from microscopies and 
presented specific results for a quasi two-dimensional XY 
symmetry magnet, providing estimates for the relevant 
length and current scales and characterizing the scaling 
of defect density as the current exceeds threshold. A 
qualitative result is that, in the presence of the dipole 
interaction, the critical current depends strongly (by a 
factor of ~ 3) on the angle between the current and 
the magnetization. Our general formalism and qual- 
itative conclusions however apply to any combination 
of dimensionality and order parameter structure which 
supports topological defects. These results complement 
the treatment of current-driven quantum criticality stud- 
ied in Refs. [HIT3I which considered transitions dominated 
by current-induced local spin excitations in Ising and 
Heisenberg magnets. 

The basic physics driving the transition, namely that 
an applied current produces a force which acts oppositely 
on vortices of different chirality, is familiar from the the- 
ory of vortices in superconductivit y 19 ' 20 and has been de- 
rive d 17 ' 18 ' 29 from the Landau-Lifshitz-Gilbert equations 
which provide a phenomenological description of mag- 
netization dynamics in the presence of applied currents. 
We note here that the relevant quantity for vortices is 
the combination of the chirality and the index giving the 
spin direction in the vortex core. 

The new contributions of the present paper include 
a derivation of the Landau-Lifshitz-Gilbert equations 
from microscopies which is based on a theory of the 
rotational degrees of freedom, complementing previous 
derivations ! 21 ' 23 ' 25 We have also included the magnetic 
dipole energy, which is found to make a numerically large 
contribution. The fact that the dipole interaction cou- 
ples spin and space directions leads to the dependence of 
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critical current on angle between current and magnetiza- 
tion. We have distinguished the transport and spin relax- 
ation scattering processes, derived an expression (follow- 
ing Ref. fill ) for the current-induced quasithermal noise, 
and explicated in detail the correspondence between su- 
perconducting and magnetic vortices. The most impor- 
tant difference is that in the magnetic case dissipative 
effects are in general expected to be weak so the vortex 
Hall angle is generically close to ir/2. 

An important question concerns the experimental ob- 
servability of the transition we discuss. Current-induced 
vortex motion has been observed^ - — The relative mag- 
nitude of the adiabatic and non-adiabatic spin transfer 
torques has been determined by studying displacement of 
vortices in permalloy discs.? 4 - In Refl32ithe nucleation and 
annihilation of magnetic vortices due to current induced 
spin transfer torque, an effect which has been theoreti- 
cally predicted ) 17 ! 18 was observed in NiFe wires. Jonietz 
et al used neutron scattering to detect a current-induced 
shift in the skyrmion lattice in the helimagnet MnSi. 
It was theoretically shown in Ref [36l and later experi- 
mentally confirmed^ that moving vortices can produce 
a "ferro"-Josephson effect where a voltage drop can be 



induced in the direction transverse to the motion of a 
vortex. The current induced vortex binding unbinding 
transition discussed in this paper should create a sim- 
ilar additional voltage drop arising from the motion of 
free vortices. Thus the effects we predict are in principle 
observable. 

A crucial simplifying assumption made in this paper 
was that the vortex-induced length scales are long enough 
that we may assume that the quasi-particle distribution 
relaxes to a steady state determined by the local magne- 
tization. It is quite likely that in many cases this locality 
assumption is violated, and that issues arise analogous to 
those arising in the physics of vortex motion in clean su- 
perconductors. Investigating this issue is an important 
direction for future research. Further important open 
questions concern the nucleation of vortices at sample 
boundaries, magnetic domain walls, and near local inho- 
mogeneities, and also the extension of our work to other 
cases such as monopoles in three dimensional magnets. 
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